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We discuss Bose-Einstein condensation of magnons (BEC) in magnets with predominant ferro- 
magnetic (FM) interaction in magnetic field H near saturation (H c ). Because H c is independent 
of FM couplings, magnetic materials of this type can have small H c that makes them promising 
candidates for experimental investigation of BEC. Ferromagnets with easy-plane anisotropy and 
antiferromagnets (AFs) containing weakly coupled FM planes or chains are discussed in detail. We 
observe small effective interaction between magnons near the QCP in such magnets, in contrast to 
AFs with strong AF coupling previously discussed. In particular, this smallness allows us to find 
crossovers in the critical temperature T C {H) oc (H c — H) 1 ^ from <f) — 3/2 to tf> = 1 in quasi-lD 
magnets, and from <j> = 3/2 to <f) « 1 (T c lnT c oc H c — H) in quasi-2D ones. 

PACS numbers: 75.10.Jm, 75.45.+J, 75.30.Kz, 75.50.Dd 

I. INTRODUCTION 

The phenomenon of Bose-Einstein condensation (BEC) is attracted much attention now. This interest is stimulated 
by recent achievements in experimental realization of BEC in ultracooled dilute gases 1 - and some antiferromagnets 
( AFs ffi & 4 -&&£ & £ ' 10 1 1 1 1 12 , It has been known for a long time that a quantum spin system is equivalent to an interacting 
Bose gasJ^ It is demonstrated in Refs i 14 : 15 that this equivalence is particularly useful for discussion of AFs in magnetic 
field H near its saturation value H c . At H > H c all spins align parallel each other and excitations of the system are 
gapped ferromagnetic (FM) magnons. As the field becomes lower than H c the gap closes and long-range AF ordering 
appears in the plane perpendicular to the field that corresponds to condensation of magnons with momentum equal to 
AF vector. The density of condensed particles is proportional to the difference H c — H that plays the role of chemical 
potential. When H c — H <C H c one can use the well-known results for dilute Bose gas. In particular, the ground 
state energy and the static susceptibility were found to differ significantly from their counterparts obtained within 
first orders of spin- wave expansion. 1 - The origin of these discrepancies is that interaction between magnons becomes 
important at H ~ H c . It is properly taken into account in Bose-gas formalism within all orders of perturbation theory. 
Consideration of the entire perturbation series is very difficult in conventional spin- wave formalism that makes Bose- 
gas description to be more efficient for discussion of this problem. Then the point H = H c is a quantum critical point 
(QCP) that belongs to BEC universality class. When H < H c the condensate disappears at a finite temperatur o 14 : 16 

T c = C(H C - H) 1 ^, (1) 

where <f> = 3/2 in 3D systems and C is a constant. 

As the saturation field in AF is proportional to the value of the exchange coupling, H c appears to be in a hardly 
achievable range in the majority of AF materials. It was the main obstacle to carrying out the corresponding 
experiments. Only recently the first experimental observations of BEC below saturation field were reported made in 
CS2CUCI4 (H c 8.5 T) ) 11 : 12 where, in particular, <j) ps 1.5 was observed^ 

In AF with singlet ground state a QCP of BEC universality class can exist in magnetic field lower than its saturation 
value^i 4 - In this case the field larger than a gap A separating singlet ground state and first triplet levels generates a 
gas of spin-triplet states moving in a non-magnetic background. These triplets can be regarded as bosonic particles. 
The difference H — A plays the role of the chemical potential that controlls the number of particles. QCPs in a 
number of AFs of this type have been investigated recently both experimentally and theoretically^^i^&L&iLiLii 

The possibility to produce BEC in quantum spin magnets has stimulated considerable experimental efforts to 
find candidate materials that have been focused on AF compounds. Then, the aim of the present paper is to draw 
attention to magnets with predominant FM interaction, which are equivalent near QCP to dilute Bose gas with pair 
interaction. Because H c is independent of FM couplings, magnetic materials of this type can have small H c that 
makes them promising candidates for experimental investigation of BEC. In contrast to AFs with strong AF coupling 
previously discussed, we observe small effective interaction between magnons in such magnets that leads to a number 
of particular features. 

We discuss in detail two classes of such magnets. One of them contains 3D, quasi-2D and quasi-lD ferromagnets 
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with an easy-plane anisotropy (EPF) which Hamiltonian has the form 

n = -\ E J ^ S i + \ E A v S * S ! + H J2 S i> ( 2 ) 

where the first term is an isotropic exchange interaction, the second one describes anisotropic exchange of the easy 
plane type with z to be a hard axis and the last term is the Zeeman energy in transverse magnetic field H directed 
oppositely to z-axis. J and A are implied to be positive and short-ranged. We assume below that planes in quasi-2D 
magnets are stacked along z axis and chains in quasi-lD magnets are parallel to z axis. H c is determined by the 
anisotropy in EPFs that is typically small being of relativistic nature. Specific calculations for EPFs similar to those 
for AF o 14 ' 15 have not been performed yet. Moreover the fact that QCP in EPFs belongs to BEC universality class 
seems to be not well recognized. To the best of our knowledge, Refs . 17 ' 18 are the only papers in which EPF near QCP 
was discussed. It was obtained for 5 = 1/2 using one-loop RG approximation in the leading order of e — 2 — d, where 
d is EPF dimension, that H c — H ~ ?c/|e| + 0(e 2 ) for d > 2, i.e., 4> = 1 in Eq. |T]). Meantime it remains unclear 
whether this formula works for |e| ~ 1 and in particular for e = —1 (3D EPF). The paper— should be mentioned, 
where XY model was discussed using RG analysis which Hamiltonian is a special case of with Ay = Jy . It was 
established there that for 2 < d < 4 the line of phase transition in H — T plane near QCP is given by Eq. (JTJ) with 
<j> = d/2. This result is consistent with our finding that the QCP belongs to BEC universality class at d = 3. 

Another class of magnets to be considered in the present paper in detail includes quasi-2D (quasi-lD) AFs containing 
FM planes (chains) weakly coupled antiferromagnetically. The Hamiltonian of such systems is given by Eq. @, where 
Jij is understood to be antiferromagnetic (negative) if i-th and j-th sites belong to different planes (chains). The 
specific calculations for this type of AFs have not been performed yet. H c is expected to be small being determined 
by small AF coupling between planes (chains) and easy-plane anisotropy Ay (if any). 

The rest of this paper is organized as follows. The Hamiltonian transformation and technique are discussed in 
Sec. HU Effective interaction between magnons near QCP is considered in Sec. IIIII Dimensional crossover in T C (H) in 
quasi-low-dimensional magnets is discussed in Sec. IIV1 Section IVl contains our conclusions. There are two appendixes 
with details of calculations. 



II. HAMILTONIAN TRANSFORMATION AND TECHNIQUE 

Let us consider first a magnet with S = 1/2. Following Refi 1 ^ we express spin projections as 

St = -\ + a\ ai , S\ = a\, Sr= ai . (3) 

where a\ and a$ are Pauli operators that commute on different sites, aj<zj + a|<Zj = 1 and (aj) 2 = af = 0. Hamiltonian 
@ is written as 

U = 22 (e k - fi) a^Ok 

k 

+ E(^ k i+ k 3 ~ J ki+k3)a kl a k 2 a - k 3 a -k 4I (4) 

where £k = 5(Jk — Jk), ko is AF vector that is equal to (0,0, it) and (w, tt, 0) in quasi-2D and quasi-lD AFs, 
respectively, and ko = in EPFs. The constant \i = H c — H plays the role of chemical potential, where H c is the 
classical saturation field. H c = SA in EPFs and H c — SA + 25| Jq| in AFs. The momentum conservation law 
ki + k2 + k3 + k4 = is implied in the sum of the second term in Eq. {I} and in similar sums below. The bare 
spectrum ek is quadratic in EPFs near k = 0: = Dk 2 in 3D EPF and £k = Dk 2 + D'k 2 ^ in quasi-2D (quasi-lD) 

EPFs, where kx is the projection of k perpendicular to the planes (chains) and k 2 = k 2 + k\. In quasi-2D and 
quasi-lD AFs the spectrum is quadratic near ko: £k+k — Dk?, + D'k\. To treat operators a as Bose ones we should 
introduce to the Hamiltonian ((4]) the constraint term U/N ^2 aj Ci aj C2 a_k 3 a-k4, where U — > oo, that describes infinite 
repulsion of particles on the same sitet^ 4 - As a result we lead to a gas of Bose particles with pair interaction only and 
with quadratic spectrum which properties are well established in the small density limit. 16 ' 20 In order to find physical 
observables one has to derive effective interaction between two particles. The ladder approximation (LA) is valid in 
3D dilute Bose gas and one has to work out the equation on the vertex T(p, 0) shown in Fig. [lja). It can be solved 
exactly as in Ref.— . Before discussion of the result let us consider larger spins. 
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An appropriate spin transformation for 5 > 1/2 was proposed in Ref. 15 basing on Holstein-Primakoff one: 

5? = -S + a\ai, (5a) 
Sj = V2Sa\ (l - (1 - K) ajfli) , (5b) 

Sr = V2S (l - (1 - K) a\a^ a u (5c) 

where K = — 1/(25). Expressions for S\ and S~ were derived from the conventional Holstein-Primakoff ones 
Sj = V2Sa\ [1 — a|aj/(25)] 1 / 2 and S~ = y/2S[l — a\ai/(2S)] 1 ^ 2 a i by expanding the square roots, putting all operators 
a\ to the left of all Oj using commutation relations and discarding terms containing more than three operators a\ 
and <Xi (it is reasonable because one can neglect interaction of more than two particles in dilute gas). In terms 
of this representation Hamiltonian @ has the form (fj| with S > 1/2 and with the additional interacting term 
5(1 — ty/NJ^iJk! + ^k3)«L«k 2 a-k 3 a_k4- Equation for the vertex T(p,0) shown in Fig.QJa) can be solved exactly 
as in Ref. 15 .. Notice that in the case of 5 = 1/2 the results for the vertex obtained using the exact method discussed 
above and using the approximate spin representations ([5]) coincide, as in AFs considered in Ref^. 



III. EFFECTIVE INTERACTION BETWEEN MAGNONS NEAR QCP 

Let us turn to 3D EPF first. The expression for the effective interaction A that is valid for 5 > 1/2 has the form 

\ EPF ee 2r(0, 0) = - A ° (6) 

J- + -lEPF/{^b) 

where Tepf = l/A r X)k^k/(^o — Jk)- Notice that A — ^ if Ay — * that should be contrasted to AFs with 
strong AF coupling, where A is proportional to the exchange coupling constant and T ~ L 14 i 15 The density of 
condensed particles po can be found minimizing the energy shown schematically in Fig. [ljb) that has the formiS 
E = -Wo + Xpl/2 + 2\p /NJ2 k N(£k), where JV{uj) = (e w / T - l)" 1 . As a result one obtains 1 ^ 

po = f-|E^)- (?) 

k 

Temperature fluctuations decrease po that disappears at a critical temperature T C (H). After simple integration in 
Eq. © we have for T c Eq. |T]) with 4> = 3/2 and C = 47rD(2C(3/2)A)~ 2/3 . 16 One finds for renormalized spectrum 
?k = [fk(ek + 2A)] 1 / 2 if A > (phase with condensate, p ^= 0) and ?k = ek + |A| if A < (phase without 
condensate, po = 0), where A = \i — 2A/AJ^ k A^(ek)J^ The average spin projections have the form (5 2 ) = —5 + p, 
(S x ) — ^pocosip, (S v ) = y/posin'i/), where tp is an arbitrary phase and p = po + 1/N A/"(ek) is the full density 
of particles. Recall that all expressions for the vertex and self-energy parts are obtained within LA discarding some 
thermal corrections J£ Meantime the results obtained are self-consistent except for very narrow area along the line of 
phase transitions, \AT\/T C <C A^/TjT/D 3 / 2 (here and below AT = T c ~ T), where the omitted temperature corrections 
become large and should be taken into account^ (see also Appendix [A] for discussions) . 

One obtains for the ground state energy E = —p 2 /(2\) while its classical value is given by E c i — — p? /{2A Q ). It is 
seen from Eq. ([6]) that A < A Q and quantum fluctuations lower the energy of the ground state as in AFs previously 
discussed i 14 ! 15 The average angle between spins and magnetic field at T = is given by 9 ~ [2p/(SX)) 1 ^ 2 while its 
classical value is 9 c i ~ [2/i/(5Ao)] 1 ^ 2 . It is seen from Eq. ([6]) that quantum fluctuations renormalize strongly the 
ground state energy and the magnetic susceptibility if T/(25) ~ 1 but one could expect T -c 1 in 3D EPFs because 
T ~ Aj J . We see that in contrast to AFs with strong AF coupling effective interaction between magnons near QCP 
is small in EPF that leads to weak renormalization of the classical values of physical quantities. One infers that this 
interaction is small also in other 3D magnets with predominant FM coupling. 

We cannot point out a compound that is a good 3D EPF without or sufficiently small anisotropy in the easy 
plane. In contrast, quasi-2D and quasi-lD materials are well-known containing FM planes or chains weakly coupling 
ferromagnetically (EPFs) or antifcrromagnctically (AFs). We obtain for the effective interaction Xaf in the last class 
of compounds 
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where T AF = 1/iV X) k (A k + 2J k )/(J ko - J k ). Eq. jSJ) is valid for 5 > 1/2, as is Eq. It is seen that in contrast 
to EPFs the effective interaction is determined by both easy-plane anisotropy (if any) and inter-plane (-chain) AF 
interaction. 

It should be pointed out that the smallness of the effective interaction A -C D in the magnets discussed requires 
to find terms proportional to to and e q in T(ko,q) at small q — (w,q). Notice that these terms are negligible in AFs 
studied in Refs . 14 ' 15 because A ~ D there. In contrast, when A <C D and T ^> A, such terms might give significant 
contribution for example to the last term in the expression for E shown in Fig. [T^b). Equation for T(p, q) presented 
in Fig. HJa) can be solved exactly in quite a straightforward manner like that for 0) (see Appendix IB1 for some 
details), the result being 



(l-£) 2 

i + r/(25)' 



2r(fco,<?) = A+(e q -o;-2 M ) T A. 7F77 ^ . (9) 



The last term here contains the difference e q — to and in consequence its contribution to the last term in Fig. [2(b) is 
negligible when 

It should be stressed that despite the formal condition of applicability of the dilute Bose-gas formalism, that is read 

1 /3 

in 3D EPF as p X/D <C 1, is fulfilled for all values of H, our consideration is valid for /i <C A only. Really, the above 
formulas lead to nonsenses at p ~ A. It is seen, in particular, that at S = 1/2 Eq. © gives r(fco, 0) s» H that results 
in unrestrictedly large density of condensed particles po — p,/(2T(ko, 0)) at small H . To explain this discrepancy we 
remind that dilute Bose-gas approach is applicable only when one can neglect effective interaction between more than 
two particles. Then, let us consider correction to the energy from diagrams containing six wavy lines that are shown 
in Fig.[ljc), where each square denotes the vertex obtained within LA. Essentially, one should take into account the 
dependence of the vertex on the external momenta to calculate such diagrams at p ~ A because summation over large 
momenta is important in this case. It is easy to realize taking into account Eq. (0) and the equality 

2r(p,0) = A + 2e p (l-/C) (10) 

(see Appendix [B] for some details of calculation of Eq. (fTTJ)0 that contributions from these diagrams are of the order 
of DpQ. Their sum is proportional to PqT^ 3 \ where is the three-particles vertex with external momenta equal 
to fco^ As a result one infers that intermediate processes of scattering of two particles at large momenta gives the 
greatest contribution to the three-particles vertex. Then we could expect that effective many-particles interaction is 

1 /3 

negligible only if p X e ff/D -C 1, where A e // ~ D is an effective two-particles interaction at large momenta. Thus, 
we lead to the same condition of applicability of the dilute Bose-gas formalism as in AFs with strong AF coupling: 

1 /3 

p <C 1. Notice that one should use effective two-particles interaction at small momenta, i.e., A given by Eqs. ^ 
and ||5J|, to find physical observables within LA at small T and p. 

The formal conditions of applicability of the dilute Bose-gas formalism in quasi-2D and quasi- ID magnets are 
(poD/D'Y^X/D < 1 and (p D 2 /D /2 )^ 3 X/D < 1 which lead at A - D to p < D' and p < D' 2 /D, respectively. At 
such small p one has T c <C D 1 and T c <C D' 2 /D, respectively, and magnets behave like 3D one with = 3/2. LA is 
valid at larger p when A -c D. It allows us below, in particular, to observe crossovers in T C (H) at p -C min{A, D'}. 
Notice that the above conditions are satisfied at such p, although it is also required for this in quasi-lD magnet that 
A < ^fDLy. 



IV. QUASI-LOW-DIMENSIONAL MAGNETS 

We demonstrate now that the observed smallness of effective interaction between magnons (A<fl and T/(2S) -C 1) 
leads to the validity of LA even at T ^S> D' in low-dimensional magnets. Notice that LA is not valid in AFs with 
strong AF exchange at T 3> D'. For some details about the estimations of the range of LA validity presented below 
(i.e., estimations of \AT\/T C ) see Appendix [A"l 



A. Quasi-two-dimensional magnets 

A magnet is considered to be quasi-two-dimensional one if \u(D/D') 3> 1. This criterion is fulfilled in such EPFs 
as K 2 CuF 4 (S = 1/2, T c = 6.25 K, J = 20 K, A = 0.8 K and J' = 1.2 • 10~ 2 K) and stage-2 NiCl 2 graphite interlayer 
compound (S = 1, T c = 18.7 K, J = 20 K, A = 0.16 K and J' = 10~ 3 K)M- KCrSe 2 is a magnet with weak AF 
coupling between FM planes (S = 3/2, T N = 40 K, J = 16.7 K, and J' » -0.06 K)^ 

When p -C XD' /D -^i H c (region I in Fig. [2j) a quasi-2D magnet behaves like 3D one and the results obtained within 
LA are valid when \ AT\/T C 3> \jT c jD'XjD. In contrast to 3D EPF, summation over small momenta is important in 
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the expressions for T and we have T EPF = SA /(AttD) ln(D/D') and T AF = S{A + 2| J^\)/{AttD) \n{D/D'). T/(2S) 
can be of the order of unity only at exponentially small D' that is not the case in the materials mentioned above. 
After simple integration in Eq. @ we lead to Eq. (J), where 4> = 3/2 and C = 4tt(.D'L> 2 ) 1 / 3 (2C(3/2)A)~ 2 / 3 . 

Properties of a quasi-2D magnet are more specific at H c ^> /i ^> XD' / D and TcfJV) ^> T ^> D'. Assuming that 
T / (2S) 1, as in all materials mentioned above, one can replace r(fco, q) by the constant A in expression for E shown 
in Fig.QTb) and use Eq. ([7|) to find T c , the result being 

T c ]n(^)«27rZ>£ (11) 

Let us discuss the range of validity of LA in this case. As in 3D Bose gas, we omit thermal corrections to the vertex 
during the calculations. Comprehensive analysis of diagrams similar to that carried out in Ref.— shows that these 
thermal corrections are negligible at H > \x > XD'/D (region II in Fig. [2]) when \AT\/T C > T c X/[DD'ln(T c /D')], 
where H ~ mm{H c , D'}. Notice that T c <C DD'/X in region II. Regions I and II exist in all compounds discussed. If 
H c 3> D' LA is valid also at H c [i ^> D 1 (region III in Fig. but the range of validity of LA is much narrower: 

|AT|/T c »l/In(T c /£»'). 

Region III is absent only in KCrSe2 among materials listed above. Unfortunately D' is too small in all these 
compounds and it seems to be difficult to reach the region I experimentally Meantime Eq. (fTT|) can be verified in 
them at T <C Term an d M <S H c . 

B. Quasi-one-dimensional magnets 



We consider a magnet to be quasi-lD one if y/D/D' ^S> 1. This condition is fulfilled in EPF (C6HnNH3)CuBr3 
(CHAB) (S = 1/2, T c = 1.5 K, J = 110 K, A = 5.5 K and J' ~ 0.1 K) 23 and in AF CsNiF 3 (S = 1, T N = 2.8 K, 
J = 23.6 K, A = 9.0 K and J' ~ -J/50Q)2i. 

Quasi-lD magnet behaves like 3D one at T <SC £)' and -C X\J D' / D <C JJ C (region I in Fig. [3]). LA is valid when 
\AT\/T C > yjT c /DX/D'. We have for T c Eq. © with = 3/2 and C* = 4tt(L>' 2 L>) 1 / 3 (2C(3/2)A)- 2 / 3 . One estimates 
Tepf ^ 5A /(47rv / iJD 7 ) and T^p ~ ^(Ao + 2|J^|)/(4 7 rVllD 7 ) and leads to A ~ A /1.5 in CHAB and A ~ A /l.7 in 
CsNiF 3 . 

Smallness of A and T /(2S) allows us to use LA at H > fi > Xy/D'/D (region II in Fig. [3|) and T C(JV ) > T > £»'. 

If T /(25) < 1 (i.e. if A < V-D£>' in EPFs and A + 2\ J' \ < V-D-D' in AFs) one can replace r(fc , t?) by the constant 
A in Fig. []Tb) and lead from Eq. to 

T c = ttVUd 7 ^-. (12) 



Analysis of diagrams shows that the range of validity of LA is given by \AT\/T C 3> T C X/ '{D 1 V DD'). Notice that 
T c < D'VDD'/X in region II. Interestingly, the condition T /(2S) < 1 holds in a quasi-lD AFs with A < \ Jq|, where 
H ~ _ff c and LA is valid at /i <Si H c as in 3D magnets. 

LA is not valid in the region II in the above mentioned compounds because T/(2S) ~ 1 in them. The exponent 
= 3/2 can be observed in these materials at sufficiently small temperatures T <C 0.1 K and at fi <§C 0.5 T. 

V. CONCLUSION 

In summary, we discuss BEC of magnons in magnets with predominant FM interaction in magnetic field H near 
saturation (H c ). EPFs and AFs containing weakly coupled FM planes or chains are discussed in detail. We find 
effective interactions between two particles in EPF © and in AF ((SJ) . In contrast to AFs with strong AF coupling ) 14 i 15 
we observe small interaction between magnons near the QCP. In particular, this smallness allows us to find crossovers 
in the critical temperature (p} from <j) = 3/2 to (/> = 1 (Eq. (| 1 2 [) ) in quasi-lD magnets, and from cj> = 3/2 to <p w 1 
(Eq. (JTTJ)) in quasi-2D ones (see Figs. [Hand E]). 
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APPENDIX A: RANGE OF VALIDITY OF THE LADDER APPROXIMATION 



We turn in this appendix to discussion of the range of validity of LA that is used for calculation in this paper. 
Firstly, let us consider phase without condensate. One has only normal Green's function G(p) = l/(u> — e^+fj, — 
where is normal self-energy part. X(p) is determined within LA by the diagram shown in Fig. 0] that gives 

2A/AT^ k 7V(ek), where we use the bare Green's function G^°'(q) = l/(w — e q + fj). As a result one has within LA 
G«(p) = i/( w - e k - |A«|), where 



2A 



(Al) 



Recall that A = gives the line of phase transitions within LA. To find in the second approximation (i.e., to 
obtain A^ 2 )) one should consider the diagram presented in Fig. [4] and use Green's function obtained within LA. As a 
result one has 



A (2) 




3D, 

quasi-2D, 
quasi- ID, 



(A2) 



where we assume that \A^\ <C D' . Expression (|A2|) should be much smaller than A' 1 ' in order LA to be valid. 
For 3D systems we can represent A near the line of phase transition as follows: |A| ~ (AT/T C )X(T C /D) 3 / 2 , where 
AT = T-T C . As a result we have from Eq. that LA is valid if AT/T C < X^/T^/D 3 / 2 . For quasi-low-dimensional 
systems corresponding estimations depend on the distance to the QCP (see Figs. [2]and[3|). In quasi-2D magnets 
one obtains the following results. Region I: |A| - (AT/T c )XTc /2 /(Dy/TJ 7 ) and AT/T C > y/T c /D'\/D. Region II: 
|A| - {AT/T C )X{T C /D) \n{T c /D') and AT/T C > T C X/[DD' ln(T c /D')]. In region III one should assume that |A| > D 1 
that gives instead Eq. (JM} |A( 2 > - ~ XT / D and we have as a result AT/T C > l/ln(T c /D'). One obtains 

for quasi-lD system in region I |A| - (AT/T C )AT C 3/2 /(D'y/D) and AT/T C > y/T c /D\/D'. In region II we have 
|A| - (AT/T c )XT c /VDty and AT/T C > T C X/[D' VtJTJ 7 ]. 

Evaluation of the range of validity of LA in phase with condensate is somewhat more tedious task because, in 
particular, one has two Green's functions, normal and anomalous, at T < T c . Meantime it can be done as in Ref^ 
for 3D dilute Bose gas. The result is that LA works below the line of phase transitions as close to the line as it does 
above the line. 



APPENDIX B: CALCULATION OF THE VERTEX 



We present in this appendix some details of calculations of the vertex. For definiteness we discuss here EPF with 
5 > 1/2, i.e., we use spin representation J5]). All calculations for 5 = 1/2 using spin representation (|3|) can be carried 
out similarly and the results coincide with those for 5 = 1/2 obtained using Eqs. ([5]). AFs with predominant FM 
interaction can be considered in a similar way (in fact, in this case one should replace Fourier components by 
Ak — 2J k in all formulas below). It should be noted also that the method of calculation resembles that has been 
proposed in Refs. 14 ' 15 to find T(p, 0) in AFs. 

The equation for the vertex in EPF with 5 > 1/2 that is represented graphically in Fig. HJa) has the form 



4T(p, q) = - Jp+ q - J P + A p+q + A p - 25(/C - 1)(J + J q + J p + J p+q ) 



-y: 



T(k,q) 



— [-^k-p + '^k+p+q — ^k-p — ^k+p+q + 25(/C — l)(Jk + Jp 



■h 



k+q 



jpim) 



where q — (w,q), p — (0,p) and we use the bare Green's function G^ \q) = l/(o> — e q + //). Let us introduce a new 
quantity 



(B2) 
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Then, we have after summation over p of both parts of Eq. (|B1[) 

4T(q) = 25(1 - K)(J + J q ) + iSQC - 1)1 £ F(M) -(J k + J k+q ). (B3) 

iv ^ e k + e k+q - lu — Zfi 

It is convenient to seek the solution for the vertex in the form 

T{p,q)=T(q)+a(q)J p , (B4) 

where a(q) is a function independent of p. Substituting Eq. (|B4[) to Eqs. (|B1[) and (|B3|) we lead to a set of two linear 
equations on T(q) and that can be readily solved. As a result one obtains, in particular, Eq. © for T(ko,q) and 
Eq. |10D for r(p,Q). 
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FIG. 1: (a) Equation for the vertex T(p, g) in the ladder approximation. Here fco = (w = 0,ko) and ko is AF vector (ko = 
in EPF), (b) Part of the energy of the dilute Bose gas depending on po- Each wavy line denotes condensed particle and 
corresponds to the factor yfpo- (c) Correction to the energy from diagrams containing six wavy lines. Each square stands for 
the vertex obtained within the ladder approximation. This correction is equal to PoT i where r* 3 ' is the three-particles vertex 
with external momenta equal to feo- 




H -H 



H-Z 



D' H H 



FIG. 2: (color online). Sketch of H — T plane for quasi-2D magnet. Three regions are marked off: n — H c — H <<£ AD' /D 
(region I), H > fj, > XD'/D (region II) and H c > /j, > H (region III), where H ~ D', f ~ DD'/X if H c > D', and 
H ~ _ff c , f < T c(Ar) if H c < D'. Solid (red) line is T c (if) that is given by Eq. JTJ with <p = 3/2 in region I and by Eq. flTTJ, 
cf> ~ 1, in regions II and III. Dashed lines restrict the shaded (blue) area inside that the ladder approximation is not valid: 
\AT\/T C < y/T c /D'X/D, \AT\/T C < T C X/[DD' \n(T c /D')] and \AT\/T C < l/ln(T c /D') in regions I, II and III, respectively, 
where AT = T c - T. 
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FIG. 3: (color online). Sketch of H 
D'VDW/X if H c > D', and H ~ H c 



— T plane for quasi-lD magnet (similar to Fig. [2] for quasi-2D), where H ~ D' , T ~ 
T Tc(n) if < -D'. The ladder approximation is not valid inside the shaded (blue) 



area determined by \AT\/T C < ^Tc/DX/D' and |AT|/T C < T C X/(D'VDD') in regions I and II, respectively. The results 
for region II are valid only if Ao <S V DD' in EPFs and Ao + 2j Jg\ <C V DD' in AFs signifying the smallness of magnons 
interaction. 
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FIG. 4: A diagram for the normal self-energy part E(p) to be taken into account in the phase without condensate. Here the 
square denotes the vertex obtained within the ladder approximation. 



